function f = myfun(x)
f = 3*x(1)^2 + 2*x(1)*x(2) + x(2)^2;


x0 = [1,1];
[x,fval] = fminunc(@myfun,x0);
